Beating in electronic transport through quantum dot based devices 
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Electronic transport through a two-level system driven by external electric field and coupled 
to (magnetic or non-magnetic) electron reservoirs is considered theoretically. The basic transport 
characteristics such as charge and spin current and tunnel magnetoresistance (TMR) are calculated 
in the weak coupling approximation by the use of rate equation connected with Green function 
formalism and slave-boson approach. The time dependent phenomenon is considered in the gradient 
expansion approximation. The results show that coherent beating pattern can be observed both in 
current and TMR. The proposed system consisting of two quantum dots attached to external leads, 
in which the dots' levels can be tuned independently, can be realized experimentally to test this 
well known physical phenomenon. Finally, we also indicate possible practical applications of such 
device. 

PACS numbers: 72.25.-b,73.23.-b,73.63.Kv,85.35.Be 
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I. INTRODUCTION 

Beating is a well known phenomenon in physics. "'^ It 
occurs when the difference between frequencies of two 
interfering waves is small enough. As a result a long- 
wavelength pattern appears (with characteristic envelope 
changing very slowly). The resulting beating frequency 
is equal to one half of the difference in original wave fre- 
quencies. This effect is very important from both the fun- 
damental and application points of view, as it provides 
a sensitive method for measuring the frequency differ- 
ence. In music, for instance, the beating effect is used 
for tuning the instruments. This phenomenon is also uti- 
lized in conventional electronics to change the frequency 
of the input signal (in so called down conversion), which 
helps to improve sensitivity and selectivity of a receiver. 
Beating effect is also used in microwave spectroscopy.- 

Recently, it has turned out that beating phenomenon 
can be observed in different quantum systems like for in- 
stance, a single quantum dot (QD)."^ Discreteness of dot's 
energy levels arising from quantum confinement make it 
able to mimic behavior of real atom, and is thus fre- 
quently referred to as artificial atom.'* Moreover, beating 
has also been reported in a qubit coupled to a fluctuator 
being in contact with a heat bathi^i^ The beating in Rabi 
oscillations'"'^^ were noticed, when the fluctuator is close 
to resonance with the qubit and the damping is weak 
enoughi^ Coherent beating in the magnetoresistance of 
ballistic tunnel junctions were also investigated^ 

The beating phenomenon in the occupation probabil- 
ity of excited state of a qubit has been predicted for 
Josephson qubit coupled resonantly to a two-level sys- 
tem (TLS), (i.e., the qubit and TLS have equal energy 
splittings).— However, this was only true when there was 
any source of decoherence. This is also why the beating 
phenomenon has not yet been experimentally verified in 
such a system. In turn, control of electron spin coherence 
in quantum dots may be provided, for instance by circu- 
larly polarized laser pulses. Consequently, quantum dots 
may enable us to observe beating. In fact, the beating 



have already been noticed in a few experiments exploring 
time dependent Faraday rotation^iiSiH in self-assembled 
QDs systems. 

However, the beating phenomenon in electronic trans- 
port through laterally confined quantum dots systems 
is an unexplored field. Moreover, there is no experi- 
ment showing beating in transport characteristics of such 
nanoscale devices. So far, investigations were mainly 
focused on the spin-independent case, where only the 
coherent oscillations were reportedJ^ii^ Recently, Souza 
has shown that the coherent oscillations become spin de- 
pendent when Zeeman splitting of the dot's level and/or 
ferromagnetic leads are consideredJ^ In this case, the 
two spin components of the current oscillate with differ- 
ent frequencies and the beating is reported for relatively 
small splitting in the frequencies (i.e., dot's level). More- 
over, Prefetto et al. have shown that intradot spin-flip 
scattering suppresses the amplitude of the beating.*^ Re- 
cently, beating in current have been predicted due to the 
presence of Andreev bound states in dc biased QD sys- 
tem (coupled to superconducting leads) and irradiated 
with a microwave field of appropriate frequencyJ^ More 
recently, the beating phenomenon in coherent transport 
through a microscale back-gated substrate coupled to op- 
tically gated quantum dot has been predicted when the 
Rabi frequencies approach the intrinsic Bohr frequencies 
in the dot 

Here, we propose another quantum system, where 
the beating can be observed. Especially, we consider 
two single-level quantum dots attached to ferromag- 
netic/nonmagnetic leads or to spin batteries. Experimen- 
tally it can be fabricated making use of a two-dimensional 
electron gas formed at the interface of semiconductor 
heterostructure. The system is designed in specific way 
to avoid the channel mixing effects between the dots. 
Thus, the indirect coupling between the dots is elimi- 
nated. Moreover, the direct hopping is also excluded 
and the dots can be treated as independent. Charge, 
spin current and tunnel magnetoresistance (TMR) are 
derived in the weak coupling approximation utilizing rate 
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equation associated with the Green function formalism as 
well as within the slave-boson approachJ^^iS The gradi- 
ent expansion is utilized to include time-dependent phe- 
nomenon. 

The paper is organized as follows. In section |IT] we de- 
scribe the model and theoretical formalism. Numerical 
results on current and TMR are presented and discussed 
in section IIIII Summary and final conclusions are gath- 
ered in section Irvl 



II. MODEL AND THEORETICAL FORMALISM 

We consider two single-level quantum dots coupled to 
external electrodes (magnetic or nonmagnetic). More- 
over, nonmagnetic leads can be driven by both charge 
and spin bias voltage. The quantum dots are attached to 
the leads as shown in Fig[T] As channel mixing effects^S 
are minimized, we are allowed to introduce two indepen- 
dent transport channels, provided some additional as- 
sumptions are also made. Specifically, we also eliminate 
direct hopping between the dots (by creation of suffi- 
ciently wide and high tunnel barrier between them). The 
indirect coupling may be significantly reduced in compar- 
ison to dot-lead coupling when, for instance destructive 
interference effects take place. In real systems such pro- 
cesses are present leading to suppression of the channel 
mixing effects. As the interdot Coulomb interactions are 
at least an order of magnitude smaller than the intradot 
Coulomb interactions, we omit the former. 

Then Hamiltonian of the system is as follows: 

kacr 1=1,2 CT 

+ J2 U^nM + E E^^-^a.?- + H.C.). (1) 
2=1,2 ka i,(T 

The first term describes here the three leads in the 
non-interacting quasi-particle approximation, where a — 
SI, S2,D means two sources and one drain leads. Here, 
''kQo- (ckQo-) is the creation (annihilation) operator of an 
electron with the wave vector k and spin a in the lead a, 
whereas ekao- denotes the corresponding single-particle 
energy. The next two terms in the Hamiltonian ([T|) de- 
scribe the two quantum dots. Here, riic ~ Qia1i<y is 
the particle number operator (i = 1,2, a =t, J,), eicr{t) 
is the discrete energy level of the i-th dot (including 
time dependence of the corresponding gate voltage) , and 
Ui is the intra-dot Coulomb integral. The last term 
of Hamiltonian ([1]) describes electron tunneling between 
the leads and dots, where 1^" are the relevant tunnel- 
ing matrix elements. Coupling of the dots to exter- 
nal leads can be parameterized in terms of Tf^(e) — 
2TTj2k^i"^ia*H^ — Ekai)- Wc assumc that Tf^ is con- 
stant within the electron band, rfg.(e) = Tf^ = const for 
e e {-W/2,W/2), and rf^(e) = otherwise. Here, W 
denotes the electron band width. 




FIG. 1. Schematic picture of two dots coupled to external 
leads. Each dot is attached to its own source lead, whereas 
the drain electrode is the same for the two dots. 



As in our model the dots are independent - they do 
not interact with each other - one can decompose the 
density matrix operator of the whole system as follows 
Ptotai = Pi ® P2 and consider each subsystem (each dot 
coupled to the source and drain leads) separately. 

Furthermore, we adopt the formalism presented in 
Ref . [l8| . Specifically, we express i-th dot's operator in 
terms of Hubbard operators^! represented by four possi- 
ble electron states in each do1i2^ which satisfied the corre- 
sponding completeness relations.^® In the next step, the 
set of auxiliary operators is introduced and the dots' op- 
erators are expressed by means of these slave-boson and 
pseudofermion operators. From the definitions of the 
Dirac brackets one is able to find the commutations (and 
anticommutations) rules for new operatorsi^ Therefore, 
the Hamiltonian of the system acquires the form: 

= E £kQO-cLaCka<T + E em(t)(/iL/2o- + d\di) 

kacT i(T 

+ E Md^ + E EKacL.(el/.. + ofU,) + H.c.]. 

i ka ia 

(2) 

Here, b\ is the slave-boson operator which creates an 
empty state in ith dot, f}^ is a peudo-fermion operator 
which creates a singly occupied state with an electron 
with spin cr, whereas d\ creates doubly occupied state 
with an electron with spin a and other electron with spin 
a in ith dot. 

In the slave-particle representation the density matrix 
elements (for each subsystem) are written in the follow- 
ing way: p^^ = e\e„ p\^ = flf.a, P22 = d\di. Here, the 
statistical expectations of the density matrix elements 
(p5j„ = (/55j„) with n = 0, cr, 2) give the occupation prob- 
abilities of the given quantum dot being empty, singly 
occupied by electron with spin-cr, or doubly occupied, 
respectively. 

To derive the rate equations we start from the von 
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Neumann equation for density matrix operator: 



(3) 



where p-j = {p'm, p\^, p'^^, P22V with j = 1,2. The ob- 
tained averaged equations for density matrix elements 
can be expressed by means of dot-lead Green functions. 
Furthermore, using Langreth theorenti^, we express the 
dot-lead Green functions by means of dot's Green func- 
tions and free leads' Green functions. After utilizing gra- 
dient expansion approximation these Green functions can 
be written in the lj space in the following way: 

GL.,dA^:t) = Vo.A9l.aG^,A^,t)+g<^^G%,iujJ)], 

Gt.kcA^^t) = v:ag:i,A^,Mc.. + G%A^,Mcal 

(4) 

where the free leads lesser (<), retarded (r) and advanced 
(a) Green functions have the following form: 



gLa =i27r/"(a;)(5(w-eka^) 



r.a 7-) 



(5) 
(6) 



with /"{oj) being Fermi-Dirac function for the a lead. 
In the above equation and in further considerations we 
omit the dot's index i as the further equations for both 
QD's acquire the same form. In Eq.(l4]) the Green 
functions of the dot, in time space, are defined as: 

G.At.t') = {{QAtUim = {{eUt)fAt)m'Hm + 

\a\^{{fUt)d{t)\S{t')Uit'))). Other parts of 
vanish for t' — t, thus are omitted as we are interested 
in t' ^ t case. Furthermore, for the sake of simplicity we 
will omit the real part in g^:^^ which is justified in wide 
band limit. Combining earlier obtained rate equations 
with Eq.Q we arrive with the rate equations expressed 
in Fourier space in the following form: 



-Fr^(i-rH)G<,(c.,i)] 

+ r^(i-rH)G<_(c.,t)-r^rHG>_(..,t) 
-r^(i-rH)G<,,(a.,t)] 

=^ / duY,[r'^n^)G>,,{^,t) 

^TAl-r{^))G<,A^M (7) 

As the transition from time space to the Fourier space 
in the time dependent phenomena is not straightforward 
it is required to justified it. Therefore, we introduced 
new time variables: a mean time i = which varies 



slowly and a fast varying time difference 5t = t ~ t', 
and expressed the Green functions in these new time 
scales, i.e. G{t,t') — > Gi6t,t)M^ Expanding G{6t,t) 
in the slow variable (F) and taking the Fourier transform 
with respect to the fast variable, we arrive at the Green 
function G{uj,t) = 2„ G(")(w, with G^") being n- 
th derivative (of the G{uj,t)) with respect to the slow 
variable. Then, we retain only the first term in this 
expansion which allows us to write the lesser dot-leads 
Green function as in Eq.Q This (lowest order) gradient 
expansion is sufficient approach as we are interested in se- 
quential tunneling regimc;^^ After exploiting the above 
obtained equations, the rate equations acquire form as 
these presented in Ref.[3 when putting intradot spin- 
fiip parameter R„ to be equal to zero. However, in the 
situation considered here, the dots' Green functions de- 
pends on both uj and the mean time i. The dots' Green 
functions wc find in the weak coupling approximation, 
deriving them from corresponding equation of motion for 
the dots operators. Technically, we assumed there is no 
coupling {Vj^ = 0) and that leads are taken to be in local 
thermal equilibrium. Thus, we obtained: 

G<^^(w,i) = i2Trp„„d{uj - €^{t}) 
Geaai'^'t) = -i2TTpoo5{oj - e^(i)) 
G^,A^, t) = i27:p226[uj - (eAt) + U)] 



G^,A^, t) = -i2^p-.-a5\^ - (e.(t) + V% 



(8) 



where the time dependence is clearly emphasized. To 
derive these Green functions we used adiabatic approxi- 
mation expanding t„ {t) around the mean time i and kept 
the terms up to linear order in the slow variable, namely 
e^ir) ~ ea-{t) + ea-{T)\T-^i{T ~t). This allowed us to write 
Jj, dtiea-{T) « ea-{t)St. Then, after making Fourier trans- 
formation, Eqs. ([8]) are obtained. Finally, connecting 
Eqs.® with Eqs. ([7]) we arrive at the coupled set of dif- 
ferential equations which we solve numerically to obtain 
time dependence of the density matrix elements. 

Current flowing from a {a = SI, S2, D) lead to the jth 
dot is obtained from the standard definition: 



(9) 



where Na is an occupation number operator in a lead. 
After performing similar calculation as above, the current 
formula becomes^^: 

J'^ = '{ j '^Y.^%ri^)\G^e.A^,t) + G>,A^,i)] 

a 

+ T%[1 - r(c.)][G< ,(c.,<) + G<,,{^,i)]. (10) 

Current passing through jih. dot can be symmetrized in 
the following way: = {J^-' — jp)/2. Total current 
flowing through the system is equal to J = + ,P. We 
assume that the distance between contacts in the drain 
lead is lesser than coherence length. 
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III. NUMERICAL RESULTS 

Time-dependent phenomenon is investigated in elec- 
tronic transport through two quantum dots coupled to 
external leads as shown schematically in FigH] The dots' 
levels are driven by time-dependent gate voltages (ac 
force), whereas the dot- lead couplings are assumed to 
be constant in time. Each dot has its own gate electrode, 
thus the dots' levels can be tuned independently. More- 
over, this allows to apply AC voltages to two dots with 
distinct external frequency and driving amplitudes. It is 
worth nothing that this can not be achieved in multilevel 
single quantum dot. 

We consider the dots' levels driven by sinusoidal AC 
voltage, and thus we assume €ia{t) = £^0- -t- SiCosilit. 
Here, Q.i is frequency, whereas 5i is amplitude of the ex- 
ternal signal applied to the ith dot. In our model, the 
chemical potentials of the source and drain leads are set 
as ^si = ^J'S2 — e.V/2 and ^id — —eV/2. Here, V is 
a bias voltage applied between the source (SI, S2) and 
drain leads. Before the time dependent signals drive 
the dots' levels, the system is in deep nonequilibrium 
due to applied bias voltage. Thus, one should expect 
the dynamics of the system undergoes non-Marcovian 
processes. In numerical calculations we assume that 
each dot is equally coupled to its pair of leads, namely, 
rf = rf = rf = rf = r with r being the energy 
unit. Moreover, we assume spin degenerate and equal 
time-independent parts of the dot levels, tia = cq (for 
i = 1,2 and a =\,\.) and equal amplitudes of the os- 
cillating signals (^i = 82 = 5). For simplicity we also 
assume the same Coulomb parameters for the two dots, 
Ui = U2 = U. 

Approximations made during calculation of the rate 
equations and current formula (gradient expansion and 
weak coupling) constrict our model to special regimes. 
There are two regimes when this approximation is valid: 

i) when k, then must be hfl <C fcsT, F ^ fcsT, 

ii) when —W/2 < a < jisi is valid for Ml W and 
F ^ W .^"^ In our numerical calculations we choose a set 
of parameters which fulfill these limitations. Moreover, 
we assume that the system initially occupies the empty 
state Pqq =1. In our calculations we also set h = 1. 



A. Nonmagnetic leads 

At the beginning we consider quantum dots coupled to 
nonmagnetic leads and assume that dots' levels are driven 
by gate voltages with different frequencies (f^i 7^ il2) but 
equal amplitudes (Si —62)- When the external frequen- 
cies differ only a little, the beating in current are observed 
as shown in Fig 15] Total current beats with frequency be- 
ing twice the difference of the frequencies of the currents 
passing through each quantum dot. Thus, the total cur- 
rent can be decomposed as a product of two parts: one 
oscillates with the average frequency F = ^{fi + f2) and 
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FIG. 2. Current as a function of time calculated for nonmag- 
netic leads and for different driving frequencies f2i = 1.0457, 
5^2 = 0.96Q. The inset shows current evolution calculated 
for nonmagnetic leads and for equal driving frequencies Qi = 
^2 = ri. Other parameters: eo = -F, S = 0.2r, (7 = F, 
fcflT = O.IF, eV = 2r. Here, Q is chosen to be frequency 
unit. 



second changing with the frequency A/ = ^{fi — f2)i 
where /i — Qi/2t: and /2 = f22/27r are corresponding 
frequencies of the currents flowing through each dot. The 
latter term controls the amplitude of the envelope and is 
responsible for the sensing of beating. The beating fre- 
quency is twice the difference frequency Fb = 2A/. Thus, 
the beating frequency is lowered when reducing the dif- 
ference in the frequencies of the input signals. This ef- 
fect is only due to the difference in the frequencies of the 
external gate voltages. To show this we calculated the 
current evolution for equal external frequencies fii = 1^2 
and displayed it in inset of Figl2l Thus, we believe that 
this system is favorable for observing current's beating in 
experiment. In contrast to results presented in Ref.|14|. 
where the beating signal is damped (due to the dot-lead 
coupling) , in our case beating of the current is sustained 
in time. 

In Fig|3]we show the influence of the temperature on 
the current's beating. We notice that the amplitude of 
the beating signal is damped as temperature increases. 
However, even for fcsT 3> F the beating pattern can still 
exist what is clearly shown in the zoomed part of Fig|31 
In turn, the amplitude of the beating can be increased by 
enlarging the amplitudes of the input signals (S). This 
implies that even for ksT > T the current's beating sur- 
vives and may be observed when S is sufhciently large. 
We also noticed that average current drops with increas- 
ing temperature, which is due to thermal damping effect 
in the leads. 

Our calculations have also shown that intradot 
Coulomb interactions do not destroy beating pattern in 
current. To show this we plot in Fig|4] beating current 
for different values of the Hubbard parameter U. How- 
ever, the Coulomb repulsion influences both the ampli- 
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FIG. 3. Beating current calculated for indicated values of the 
temperature. Other parameters as in Fig[2l In the ellipse we 
show zoomed part of the plot encircled in the frame. 



tude of the beating and the value of the average current. 
Namely, when there is no Coulomb interaction, the am- 
plitude of the beating is most pronounced. When the on- 
dot Coulomb repulsion is present, the amplitude of the 
beating is suppressed. The dependence of the beating's 
amplitude is nonmonotonic function of the parameter U. 
When Co + U is within the transport window it decays 
with increasing J7, but for eo + U > starts to in- 

crease. However, it never again reaches the maximum 
value. 

The average current is a also nonmonotonic function of 
the Coulomb parameter U. It reaches high values when 
Co ~ M-D and eg + U approaches fisj (but not very close, 
eo + U ^ ^Sj, due to imposed gradient's expansion con- 
dition.) When eo + U is beyond the transport window, 
average current drops and saturates for sufficiently large 
U . Moreover, Coulomb interactions introduce small hor- 
izontal asymmetry in the beating pattern. 
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50 100 150 200 250 

nt/27r 

FIG. 4. Beating current calculated for indicated values of the 
intradot Coulomb repulsion parameter. Other parameters as 
in Fig[2l Here, for clarity we rescaled the current plot for 
U = T and U = 4r. To obtain calculated values of the current 
for those parameter one has to add certain value to all points 
as is indicated by arrows. 




FIG. 5. Magnetic configurations taken into account. 



B. Ferromagnetic leads 

When, the leads are ferromagnetic several magnetic 
configurations are possible. To 'measure' the difference 
in these distinct configurations it is convenient to in- 
troduce tunnel magnetoresistance (TMR). This quan- 
tity results from spin-dependent dot-lead tunneling pro- 
cesses, which in turn leads to the dependence of trans- 
port characteristics on magnetic configuration of the sys- 
tem. The TMR is quantitatively described by the ra- 
tio TMR = Ip — Iap/Iap^ where Ip and Iap denote 
the currents flowing through the system in the paral- 
lel and antiparallel magnetic configurations, respectively. 
Introducing the spin polarization pa of lead a {a — 
Sl,S2,D) as Pa = {pt - Pa)/{pt + Pa), the coupling 
parameters can be expressed as T^^ = Ti^i^ ^ Pa), 
with Via = (r+ -I- r~^)/2. Here, p+ and p~ are the den- 



sities of states at the Fermi level for spin-majority and 
spin- minority electrons in the lead a, while and F~j 
describe coupling of the ith dot to the lead a in the spin- 
majority and spin-minority channels, respectively. 

Let us first consider the case where the magnetic mo- 
ments of the source leads are pinned (with 'up' direc- 
tion) and the magnetization of the drain electrode can 
be changed from 'up' to 'down' as schematically is shown 
in Figinja). We calculated the currents in both mag- 
netic configurations and TMR for leads' polarization 
Psi = PS2 = Pd = 0.5. Firstly, one observes that the 
beating is still present in the current characteristics for 
both magnetic configurations. However, TMR exhibits 
beating pattern a little distorted. As the dots are de- 
coupled from each other and we consider the case of 
weak couplings we should expect positive TMR, which 
is clearly displayed in FiglHKb). Then, off course, the 
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FIG. 6. a) Current as a function of time calculated for two 
magnetic configurations from Fig[S][a) and for p — 0.5, U — 
4r. b) Time evolution of TMR. Other parameters as in Fig[2] 
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FIG. 7. Time dependence of the TMR calculated for magnetic 
configurations from Fig[S][b). Other parameters as in Fig|S] 



current in parallel magnetic configuration is greater than 
that in antiparallel one (see FiglHlJa)). However, spin 
symmetry breaking processes, as spin-flip scattering, may 
change the sign of TMR as shown in Ref. [l^ . Here, we 
do not consider such processes. Recent experiments have 
shown that the spin relaxation time in quantum dots can 



reach millisecond^^ — or even second timescales which 
is much longer than electron tunnel rate (~ r~^). 

Now, we consider the situation when the magnetiza- 
tion of the drain lead and one of the source electrode 
are pinned, whereas the magnetic moment of the sec- 
ond source lead can be flipped, as is shown in FiglSJb). 
In this case, the difference between 'parallel' and 'an- 
tiparallel' configurations is less visible, which results in 
suppression of TMR. However, the oscillating character 
is still conserved and here the beating is even more pro- 
nounced (see Figl?]). The suppression of the TMR for this 
magnetic configurations is clear when one notices that for 
this case only one transport channel is partially blocked 
(due to relevant difference in the orientations of the leads 
magnetic moments in the AP configuration), whereas in 
the former case both channels are bad 'conductors' in the 
AP configuration. Moreover, in this magnetic configura- 
tion, tt/2 phase shift is induced in the TMR pattern. A 
small distortion in the upper semicircles comes from the 
different symmetries of the current profiles for P and AP 
configurations in the vicinity of the node points. 



C. Spin-biased leads 

Here, we consider the double dots' system subjected to 
the source and drain spin batteries^"— which provides 
pure spin current without accompanying charge current. 
Pure spin current is one of the most important points 
for spintronics. However, so far spin control methods in 
commercial devices mainly have relied on usage of mag- 
netic field^^i^ or optical techniques which are not very 
efficient. Recent experiments show that pure spin current 
can be ail-electrically generated in a micron-wide chan- 
nels of a GaAs two-dimensional electron gas^^*^. This 
is very important from the application point of view, be- 
cause other quantum systems can be easily integrated 
with such all-electrically controllable spin battery. To 
control such a device we do not need optical or magnetic 
fields which precisely adjusting is rather great effort and 
thus, useless for commercial applications. 

First, we investigate DQD system connected to sym- 
metric dipolar spin batteries, i.e. we assume that fJ-sjt — 
^Di and jJ^Sji — fJ'D't for (j = 1,2). Introducing the spin 
bias Vs, generally we may write fisja- — g{V -\- aVs)/2 
and /x_dct = —e{V + <tVs)/2 with a — 1 {a — —1) for 
a =t (cr =4-)^. As we are interested in pure spin current 
we further set bias voltage equal to zero = 0. In this 
case the net charge current vanishes, because all spin- 
up electrons flow in one direction and equal amount of 
spin-down electrons flow in the opposite direction, and 
only pure spin current is generated. The spin current is 
defined in the following way, Jg — {J^ — J^) / where 
Jcr = + (ct =t,i). However, it is worth to mention 
that when the dot's energy level is split i.e., e-|- 7^ ej,, the 
nonzero charge current can be generate d'^^i'^^'^^ . In FigU] 
we show time evolution of the spin current calculated for 
different strengths of the intradot Coulomb interactions. 
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FIG. 8. Time evolution of the spin current calculated for 
indicated values of the intradot Coulomb repulsion parameter 
and for eo = in case of the symmetric spin batteries. Other 
parameters as in Fig[2] 



One can notice that the spin current exhibits more com- 
plicated beating pattern (similar as TMR in Fig|6]). On 
the other hand, for noninteracting case {U = 0) we notice 
the clear evidence of pure beats in the spin current. How- 
ever, for this case the amplitude of the beating is small. 
When Coulomb interactions are turned on, symmetric 
beating pattern vanishes and even more features appear. 
As intradot repulsion increases, the beating in the spin 
current become more and more asymmetric and the node 
points cease to exist. Instead of node new oscillations 
emerge. Namely, for nonzero U, spin current evolution 
composes from two kind of oscillations: main oscillations 
and some sub-oscillations emerged in the vicinity of the 
nodes points (existing in noninteracting case). For suffi- 
ciently large U, the main oscillations become very asym- 
metric and the sub-oscillations are more pronounced. In 
contrast to the nonmagnetic case, the dependence on the 
Coulomb repulsion is here much more complex. The am- 
plitude of the beating is small for both small and large 
enough U. This is because the state eq + U is far away 
from the chemical potentials of the leads. However, when 
eo + U is outside the transport window, the average spin 
current grows meaningly (in contrast to the charge cur- 
rent in nonmagnetic case from FigH]). When U is suffi- 
ciently large (i.e. eo + U :s> Haa) the probability of double 
occupancy drops almost to zero (p22 ~ 0) and the occu- 
pation numbers n„ also decrease (however increases) . 
This enables effectively faster tunelling processes through 
QD and thus enlarges the spin current. For U ^ the 
spin current becomes saturated. 

Now, we consider DQD system attached to asymmet- 
ric spin batteries. In this case we set ^isja = <yeVs 
and ^iDa — 0. Let us first consider noninteracting case 
([/ = 0). The dots' energy levels are situated symmetri- 
cally with respect to the spin bias voltages of the source 
leads, e.g., Cq is in the mid between p^sjt ^-^d ^^Sj^.■ Thus, 
the same amount of spin-up electrons flows in one direc- 




nt/2ii 

2 4 6 8 

U/F 

FIG. 9. Stationary charge (a) and spin (b) current as a func- 
tion of the Coulomb repulsion parameter U. Insets: Time 
evolution of the charge and spin current,respectively, calcu- 
lated for for £0=0 and for U — 0, in case of the asymmetric 
spin batteries. Other parameters as in Figl2] 



tion and equal amount of spin-down electrons flows in 
the opposite direction and the average spin current is 
nonzero. However, due to oscillations of the dots' levels, 
the charge current is also generated, but on average it 
vanishes. In the case of asymmetric spin batteries both 
the spin and charge current exhibits well-defined beat- 
ing pattern as shown in the insets of FiglHl When the 
Coulomb interactions are turned on, a nonzero average 
charge current is induced. This is because earlier men- 
tioned symmetry is now broken. It is worth noting that 
such symmetry exists also when eq — However, 
for eg < /iD and sufficiently large U the system is in 
the Coulomb blockade, thus, we expect zero current in 
the weak coupling regime. Firstly, for a small value of 
the parameter U, the average charge current grows very 
fast reaching maximum value for U ~ 0.4r and then is 
unchanged with further increase of the [/, until eo + U 
exceeds /isj-fj when it becomes reduced a little and satu- 
rates. This drop in average charge current is because the 
state eo-l-C/ ceases to contribute to the transport. In turn, 
the average spin current, generally, grows with increas- 
ing parameter U (regardless a certain ranges of U where 
the average spin current is constant). When U is suffi- 
ciently large the average spin current is also saturated. 
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To show this dependencies we plotted stationary charge 
and spin current in Fig jS] which may be regarded as aver- 
age values of respective currents in time-dependent phe- 
nomenon. However, one should bear in mind that this 
is not true for cq + U being close to fist "^^^ to gradi- 
ent expansion condition. Hence, this range should not 
be disregarded. 

It is also worth noting that for cq < /i_D both the av- 
erage charge and spin current can change the sign. As a 
result one should expect negative charge and spin differ- 
ential conductances. Moreover, in contrast to the sym- 
metric spin batteries, here, the beating structure in spin 
current is very symmetric for all values of the Coulomb 
interactions parameter U . 



IV. FINAL CONCLUSIONS 

In summary, we have studied coherent transport 
through two uncoupled quantum dots, which are at- 
tached to nonmagnetic and/or ferromagnetic leads. Gen- 
erally, two magnetic configurations were discussed. We 
took into account the Coulomb interaction between elec- 
trons on the same dot and calculated transport char- 
acteristics in the nonlinear response regime, using the 
rate equation approach connected with Green functions 
method and with slave-boson formalism. Our analysis 
was performed for oscillating dots' energy levels within 
the gradient expansion approximation. 

We have found clear evidence of both charge and spin 
current beating as well as the beating pattern in TMR. 
We have shown that the effect is due to the difference in 
the frequencies of the applied gate voltages to the two 
dots. In magnetic case, beating in spin current or TMR 
may be deformed. However, for DQD system coupled to 
the asymmetric spin battery spin current exhibits well- 
defined beating structure. 



In this paper we have omitted the interdot Coulomb 
repulsion as in real systems it is much more smaller than 
the intradot Coulomb interactions. Moreover, for the 
parameters assumed in this paper the interdot Coulomb 
integrally would be (much) lesser than the dot-lead cou- 
pling strength, and that's why it does not lead to the 
splitting in the dot's density of state. Correspondingly, 
sufficiently small interdot interaction does not affect con- 
sidered phenomenon and is irrelevant. However, suffi- 
ciently strong interdot Coulomb interaction can intro- 
duce some deviation in the beating pattern. 

The proposed system can be used as a device to mea- 
sure frequency of an unknown signal. Then, one needs 
only one QD in one arm coupled to the source and drain 
leads, whereas the second arm delivers the unknown sig- 
nal. The arm with QD plays role as the reference channel, 
and thus, tuning the frequency of the reference signal one 
is enable to detect the frequency of the 'unknown' signal. 
Moreover, the DQD device presented above can be uti- 
lized in coding information (signal). Thus, such device 
may be called nanoscale superheterodyne. Using such a 
device we are able to mix two signals of slightly different 
frequencies. As a result, one obtain resultant signal being 
a composition of slow- varying and fast- varying parts (as 
mentioned in Sec lIII A]) . Then, one of the signals with, 
for instance low frequency, may be extracted and further 
processed. The advantage of the device is that a signal 
with lower frequency is easier to be processed. 
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